Oscillating Casimir force between impurities in one-dimensional Fermi liquids 



J.N. Fuchs, 1 A. Recati, 2 and W. Zwerger 3 

1 Laboratoire de Physique des Solides, Universite Paris-Sud, CNRS, UMR 8502, F-91405 Orsay, France. 
2 CRS BEC-INFM, Povo and ECT , Villazzano, 1-38050 Trento, Italy. 
3 Technische Universitat Munchen, James Franck Strasse, D-85748 Garching, Germany. 

(Dated: February 6, 2008) 

We study the interaction of two localized impurities in a repulsive one-dimensional Fermi liquid 
via bosonization. In a previous paper [Phys. Rev. A 72, 023616 (2005)], it was shown that at 
distances much larger than the interparticle spacing the impurities interact through a Casimir-type 
force mediated by the zero sound phonons of the underlying quantum liquid. Here we extend these 
results and show that the strength and sign of this Casimir interaction depend sensitively on the 
impurities separation. These oscillations in the Casimir interaction have the same period as Friedel 
oscillations. Their maxima correspond to tunneling resonances tuned by the impurities separation. 
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Interacting one-dimensional fermions are nowadays 
available in many experimental systems. They can be 
electrons in semiconductor quantum wires, in carbon 
nanotubes, in conducting polymers, etc. [1] or fermionic 
atoms in ultra-cold gases [2] . Impurities have a profound 
effect on the physical properties of these low-dimensional 
systems. In certain exceptionally clean systems, impu- 
rities can be created and controlled. For example, it 
was recently proposed that atomic quantum dots - i.e. 
individually trapped atoms playing the role of impuri- 
ties - could be created in an ultra-cold atomic gas 
Another possibility rests in the "buckles" artificially cre- 
ated with AFM tips in single- wall carbon nanotubes and 
which behave as tunnel junctions Q. As a last example, 
we mention atom wires that were recently created on a 
Si(553)-Au surface and which featured vacancy-like de- 
fects that were manipulated with an STM tip [5||. Once 
it is recognized that impurities can be manipulated, one 
can ask new questions in which they appear as main 
actors. For example, what is the interaction between 
impurities which is mediated by the quantum liquid of 
fermions? This problem was recently addressed in 0. It 
was found that when the impurities are far apart they be- 
have like acoustic mirrors and interact via a Casimir-like 
force. This is the expected behavior for any system whose 
low energy properties are described in terms of phonons. 
However, for the special case of non-interacting fermions, 
the calculation can be done exactly and it was shown that 
this interaction features oscillations as a function of the 
impurities separation (f|. In the present paper - part 
II - which can be considered as a continuation of JO] - 
to which we shall refer as part I - we extend some of 
these results to interacting fermions. Most significantly, 
we show that the oscillation phenomenon in the interac- 
tion remains essentially unaffected in the case of interact- 
ing (repulsive) fermions. In the light of proofs showing 
the absence_of repulsive Casimir interactions for the pho- 
tonic field this is a quite remarkable situation. Here 
we have a Casimir force whose strength and sign can be 
tuned by the impurities separation. 

In order to investigate that problem, we consider the 
following model: two localized delta impurities embed- 



ded in a repulsive one-dimensional Fermi gas, whose low- 
energy behavior is that of a single channel Luttinger liq- 
uid with parameter K < 1. The interaction between 
the impurities - which is mediated by the quantum fluc- 
tuations in the liquid - is considered in the adiabatic 
limit in which changing the distance between the impu- 
rities is slow compared to the time the system takes to 
equilibrate. Within such a limit, the interaction energy 
directly follows from the groundstate energy. We should 
stress that we focus on interactions mediated by quantum 
zero-point fluctuations and not on classical effects due to 
modification of the average particle density. The latter 
dominate at short distances of the order of the inter- 
particle distance and lead to a kind of "mattress effect". 
For a discussion of the classical ground state energy, see 
part I. 

We first consider non-interacting fermions and discuss 
the oscillations in the interaction energy in that case. 
Two delta impurities located at X\ = and a; 2 = r 
are embedded in a one dimensional (ID) Fermi gas with 
Fermi velocity vf — ^Po/m where po = kp/ir is the av- 
erage density. The impurities interact with the particles 
through a potential g a S(x — x a ) where a = {1,2} and 
g a are microscopic scattering amplitudes. When the dis- 
tance between the impurities r is much larger than the 
average inter-particle distance 1/po, it was shown in part 
I that the renormalized interaction energy between the 
two impurities is given by 
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In the preceding equation, Li 2 is the di-logarithmic func- 
tion, 3? denotes the real part, 7 Q = g a /vF is the dimen- 
sionless impurity strength and we introduced the phase 
x( r , 7a) such that x = 2fcpr + S modulo 2ir and — n < 
X < 7r. The phase shift 5(j a ) = — arctan7i — arctan7 2 = 
— arctan[(7i + 7 2 )/(l — 7172)] reflects the time delay of 
fermions scattering on delta impurities. As the function 
V\%/{vf/t) is 27T-periodic in x, the interaction energy Vu 
oscillates as a function of r. Maxima are separated by 
A(2kpr) = 2tt, which corresponds to the average inter- 
particle distance Ar = 1/po, just as Friedel oscillations. 
Maxima occur when cos(x/2) = cos(kpr + 8/2) = 0. 
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Figure 1: Interaction energy V12 between two identical strong 
impurities separated by a distance r. Dashed line: non- 
interacting fermions Vrz = i>FKLi2(— e lx )/2nr. Full line: in- 
teracting fermions (with Luttinger parameter K = 0.6) V12 = 
— ivvF/24:Kr + vfX 2 /^rK 2 . In both cases, x — 2kpr + 5. For 
non-interacting fermions and strong impurities (7 3> 1) the 
phase shift S = — 2arctan7 » —tt. For interacting fermions 
(K < 1) and strong impurities (7 S> 1), we expect the phase 
shift to be S = — 2 arctan(A" 2 7) ~ — n, see the main text. The 
full line corresponds to this guess 5 ~ — n. 



This is precisely the condition for tunneling resonances 
in a non-interacting ID Fermi g clS , clS shown in Appendix 
[Al In Appendix [Bj we discuss an interesting situation in 
which the oscillations can disappear. 

In the limit of strong impurities (7^ ^> 1), the phase 
shift S w — 7r and Eq. (TJ becomes 
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where x ~ 2kpr — n. The last equation applies to 
— 7r < x < n with a periodic extension beyond that. 
The first term has exactly the form of a ID Casimir en- 
ergy — Tvvp/2Ar, with the Fermi velocity vf playing the 
role of light velocity, see e.g. @. The second term in- 
troduces oscillations. Indeed, the function Vyz/ (vp /r) 
is 27r-periodic in x an d features wide maxima made of 
pieces of parabolas and centered at x — n , see Figure [U 
On the one hand, at a minimum, where x = 0, the inter- 
action energy is given by the Casimir energy — irv F /24r. 
This correspond to a situation of no resonant tunneling 
(anti- resonances) . The Casimir energy has the usual in- 
terpretation in terms of phononic "radiation pressure" on 
the impurities: the outer radiation pressure of phonons 
overcomes the inner one and hence the resulting attrac- 
tion between the impurities. On the other hand, at a 
maximum, where X = tt, the interaction energy changes 
sign and becomes positive reaching irvp/12r. A physi- 
cal understanding of this change of sign in the Casimir 
force might be obtained by noting that in a situation of 
resonant tunneling (% = n), fermions approaching the 
double barrier structure from outside can tunnel with no 
reflection. There is thus no contribution to the radiation 
pressure due to the phonons in the region outside the 
impurities. By contrast, fermions which come from in 



between the two barriers, have to tunnel through a sin- 
gle impurity and do not benefit from resonant tunneling. 
At low energies, they are almost perfectly reflected and 
therefore contribute a positive radiation pressure on the 
impurities. At resonance therefore, the inner radiation 
pressure is larger than the outer one, and thus the effec- 
tive interaction between the impurities is repulsive. In 
a sense, this is like having impurities with transmission 
properties that are not right-left symmetric and makes a 
connection to Imry's idea of tuning the sign of the QED 
Casimir effect using asymmetric mirrors [3|. 

We now come to the main point of part II. We 
show that the basic oscillation phenomenon contained in 
Eq |f2|), and found in part I for non-interacting fermions, 
remains essentially intact when including interactions be- 
tween the fermions. In order to study ID interacting 
fermions, we use the Luttinger liquid formalism [l(J and 
concentrate on repulsive interactions. For simplicity, we 
consider spinless fermions, for which the low-energy de- 
scription is given by the following hydrodynamic action 
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where (3 = 1/T is the inverse temperature and 1 
to the density of the liquid by 

p(x) = (p + d x 6/ir) [1 + 2cos(26> + 2k F x) + 
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where pa is the equilibrium density and only the first 
harmonics are written [ltj • The Luttinger liquid descrip- 
tion depends on three parameters: the sound velocity 
u, the Luttinger interaction parameter K and a cutoff 
energy uj c ~ uk,F above which this effective description 
breaks down. Repulsive short-range interactions imply 
< K < 1. With long-range (Coulomb) interactions, the 
system becomes a ID Wigner crystal, which can be seen 
as the K — > limit of the Luttinger liquid [ll|, [I2]> but 
can not be treated within the present formalism. Non- 
interacting fermions correspond to K — 1 and u = vp. If 
there is full translational invariance uK = vp [10]. 

We consider, as before, two localized impurities sitting 
at position x\ = and X2 = r. The interaction between 
the impurities and the Luttinger liquid is 
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where g a is a phenomenological scattering amplitude. 
As the Luttinger liquid is an effective description valid 
at low energy, the phenomenological scattering ampli- 
tude 7j a is not equal to the microscopic scattering am- 
plitude g a . The precise relation between the two scat- 
tering amplitudes is not generally known. In the renor- 
malization group language, the phenomenological scat- 
tering amplitude g is roughly a running scattering am- 
plitude at the energy scale oj c . For simplicity, we con- 
sider the symmetric case ffi = 52 = 9- Inserting the 
density expansion |(4]) into the interaction §5§ gives rise 
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to several different terms. The first term 2Pgpo is just 
twice the classical energy cost of adding an impurity to 
the liquid: it does not contribute to the renormalized 
interaction energy V12. The second term due to for- 
ward scattering is proportional to d x 9. This term is 
often eliminated via a gauge transformation: 9(x) — » 
9{x) - Kg(H(x - r) + H(x) - l/2)/«, where H(x) is 
the Heaviside step function. However, this operation 
is not totally innocent here and forward scattering can 
not be eliminated: the gauge transformation gives rise 
to a phase shift S = —2K~g/u in the 2fcp backscattering 
term. In the limit of weak impurities, this phase shift 
agrees with the one found for non-interacting fermions 
5 = — 2 arctan(g/i>i?) ~ —2g/vF- Actually, we expect 
that the phase shift in the interacting case is given by 
S = — 2 &rct&n(Kg/u) for arbitrary impurity strengths, 
however to our knowledge there are no exact results on 
this problem. In fact, the precise value of S depends on 
the detailed potential and the behavior of the fermion 
liquid at high energies where a LL description no longer 
applies. This issue is often discussed in the context of 
the x-ray edge singularity or in the Kondo problem and 
is related to the order in which two non-commuting lim- 
its are taken: sending the size of the impurity to zero 
versus sending the bandwidth to infinity. When using 
the Luttinger liquid formalism, one automatically takes 
the second limit first, see Ref. The third term due 

to 2k p backward scattering gives the dominant contribu- 
tion to the action Si and, when including the phase shift 
6, appears as 

Si&V / dr(cos[20(O)] + cos[2(9(r) + X ]) , (6) 
Jo 

where V — 2gp and \ ls a phase such that x = 2kpr + 5 
modulo 2tt and — ir < X < tt. There are also higher or- 
der terms coming from the density harmonics, however 
their naive scaling dimension shows that they are less 
relevant and we therefore neglect them. Actually pertur- 
bative renormalization group calculations show that 2kp 
backscattering is always relevant as long as K < 1 while 
all higher backscattering operators cos(2n8(x) + 2nkpx), 
with n = 2,3,..., are irrelevant as long as K > 1/n 2 
[3, [HI]. As we neglect those higher backscattering oper- 
ators, our model only describes the situation of two impu- 
rities in an interacting Fermi gas such that 1/4 < K < 1 . 

As a side remark, we note that there are other ways 
of coupling impurities to a bosonic bath of oscillators, 
such as a Luttinger liquid. For example, in the context 
of QED Casimir force studies [if}], the impurities (mir- 
rors) are usually taken to couple directly to the density 
of the photonic field, i.e. 9 2 rather than cos 29 would 
appear in equation |(6|). This is the case in particular for 
the theorem of Kenneth and Klich on the absence of a 
repulsive Casimir force Q- Another example was consid- 
ered in Ref. 0], where only forward scattering was kept, 
resulting in a mapping to the spin-boson model. These 
different models, which all describe two impurities in a 
bath of quantum oscillators, do not necessarily show the 
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same behavior. 

From now on, we shall work with the action S = So+Si 
with Si given by Eq. JfT]). The partition function is Z = 
fT>6(x,T)exp(—S) and the integral on the field 9(x,t) 
can be performed everywhere except at the position of 
the impurities x = x a . In this way we are left with an 
action A = A + Ai for the fields 9 a (r) = 9(x a ,r) only 

A ° = H^(l^l,n| 2 + l^,r l | 2 -e Il (0 1 ,_„^„+C.C.)) 
f 

A, = V / dr(cos[26» 1 ] + cos[26» 2 + X }) , 

where uj n are bosonic Matsubara frequencies, the Fourier 
transform is defined as 9 n — T J dr exp(iw„r)6'(r), the 
characteristic finite-size frequency is ui r = u/r and we 
defined e„ = exp(-|w„|/w r ) and /„ = /3|w„j/(l - e£). 
Up to an overall normalization factor, the partition func- 
tion is now given by Z — J T>9 a {T) exp(— A). The 
partition function gives the free energy F = —TlnZ, 
from which the renormalized interaction energy follows: 
Viz(r) = F(r)—F (00). From the parametric dependance 
of the action A on the phase x, we see that the free en- 
ergy (and therefore V12) will be 27r-periodic in x- As 
X = 2kpr + S, this means that V12 will have some struc- 
ture with periodicity 1/po in r just as for non- interacting 
fermions. Therefore, the interaction energy V12 oscillates 
as a function of the impurities separation r. 

We first evaluate the partition function using the sad- 
dle point approximation. For small K and strong impu- 
rities 7 = g/vF 1, the fields 9 a are pinned to the val- 
ues that minimize A4, We therefore expand the cosines 
around a minimum to quadratic order in the fields. It 
should be noted that terms independent of the fields give 
an important contribution here [20]. The resulting free 
energy is: 
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As the validity of the saddle point approximation requires 
that 2-kKV ^> u) r , the renormalized interaction becomes: 
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In the limit of zero temperature uj r » T, it gives 
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where we used uj c 3> uj r [2lj . This result generalizes 
Eq. ([2]) when interactions between particles are taken into 
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account. It also admits an interpretation as an attractive 
Casimir energy —nw r /24 - albeit with the true sound 
velocity u replacing the Fermi velocity Vp - modified by 
oscillations contained in the term UrX 2 /8irK , see Figure 
[TJ The width of the maxima measured relative to their 
amplitude is proportional to K. Therefore the maxima 
are narrower in the presence of interactions between the 
particles. The exact position of the maxima is given by 
the values of r such that x = 7T, i.e. cos(x/2) = cos(fcpr+ 
5/2) = 0, and depends on both the impurity strength and 
the inter-particle interactions through the phase shift 5. 
As in the case of non-interacting fermions, this is exactly 
the condition at which tunneling resonances occur in a 
Luttinger liquid, see the Appendix lAl 

At this stage, it is worth mentioning that there is a 
mapping of our problem - as given by actions |(3]) and 
([6]) - on the double-boundary sine-Gordon model studied 
in the context of conformal field theory [TtJ ■ For special 
values of the Luttinger parameter K = l/[2(l+n)], where 
n is a strictly positive integer (the so-called reflectionless 
points), the authors of Ref. [17] were able to calculate 
exactly the groundstate energy of this model [12] • In the 
limit of strong impurities and for the reflectionless points, 
it yields the same interaction energy V12 = — 7rw r /24 + 
uJrX 2 i l + n)/An as Eq. flQ]). 

In order to study what happens for values of K close 
to one, we turn to the self-consistent harmonic approx- 
imation (SCHA), which is just an implementation of 
Feynman's variational principle, see, e.g., Expand- 
ing the cosines in Eq. to second order in the fields, 
we replace the coefficient 2V in front of the quadratic 
term by a variational parameter A. The true action 
S = Sa + Si is now replaced by a quadratic variational ac- 
tion S v . The parameter A is determined self-consistently 
by finding a minimum to the variational free energy 
F V (X) = —Tin Z v + T(S — S v ), The latter is given by 
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the second equality being justified when u> c u> r 3> T 
and nKX ^> cu r which is a validity condition of the SCHA. 
Extremizing the variational free energy gives: 



A = 2^x P (-!) = 2^-^- 
oX \irKX + Lu c 
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When K < 1, this equation has two solutions: (corre- 
sponding to a maximum of F v ) and A > (correspond- 
ing to a minimum), where A » 2V when 2irV 3> oj c and 
A « 2V(2nKV/u £ ) K ^ 1 - 1 ^ when 2irV < uj c [H]. As 
neither A nor df(X)/dX depends on r, the renormalized 
interaction energy is given by V12 ~ /(A, r) — /(A, 00) w 



— iru> r /24: + ui r x 2 /SttK . The result obtained in the saddle 
point approximation Eq. fTO)! is recovered and is there- 
fore valid more generally for K < 1 as long as irKX ^> uj r , 
i.e. for strong impurities or large distances between the 
impurities. 

The situation of weak impurities or short distances 
(but always in the limit r> 1/po) can be studied thanks 
to a perturbative calculation of the partition function. In 
the regime K > 1/2, the interaction energy is 



V12 = -7 
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where T(x) is the gamma function and \ ~ 2kpr when 
7 <S 1, as already shown in part I. This is a RKKY- 
type interaction, which decays with a renormalized power 
law r 1_2K . It is valid as long as | Vx2 j "C ui r , i.e. r <C 
Pa 1 j 1 /( K ~ 1 K In a similar way, we obtain 
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in the regime < K < 1/2 and 



V12 = -7 cosx x 



ln(/3w r ) , (15) 



for the special case K = 1/2, again provided that 
uj r 3> T. For < K < 1/2, the interaction given by 
Eq. (fl4"ll still shows Friedel oscillations but its envelope 
does not depend on the inter-impurity distance r any- 
more, which reflects the tendency of the system toward 
Wigner crystallization when K goes to zero. However 
the validity of this calculation is restricted to small but 
finite temperature as | V12 1 <C uj r corresponds to r < 

0, 



ulu^T /(vfIsfj) 1 ■ This means that when T 
the domain of validity of perturbation theory shrinks to 
zero, when < K < 1/2. We note that, at the reflection- 
less point K — 1/4, Eq. (fl4"ll has the same structure as 
the exact T = result V12 — — 7 2 cosx x (w^fc^) 2 / (iruj c ) 
where x ~ 2kpr when 7 <C 1. Again, as in the case 
of strong impurities, maxima in the interaction energy 
V12 occur when cos(x/2) = 0, which is the condition for 
tunneling resonances in a Luttinger liquid, see the Ap- 
pendix lAl 

Tunneling resonances for the conductance are only 
found for 1/4 < K < 1. Indeed on resonance 2fcp 
backscattering is absent and higher order backscatter- 
ing operators control the behavior of the system. Once 
K < 1/4, 4&f backscattering, which we neglected, be- 
comes relevant and washes out the tunneling resonances 
[l4T |. Therefore, in order to study what happens for 
very small values of K, one should also take these other 
backscattering operators into account when computing 
V12. Using the methods of part I, the present work could 
also be extended to treat spin 1/2 fermions. 

In conclusion, we found that at large impurities sepa- 
ration or for strong impurities, the renormalized inter- 
action between the impurities results from a Casimir 
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force modified by oscillations - due to quasi-stationary 
states in between the impurities. The oscillations can 
be spanned by changing the distance between the im- 
purities and are separated by the average inter-particle 
distance, i.e. they have the same periodicity as Friedel 
oscillations. The position of the maxima (peaks) is given 
by precisely the same condition as for tunneling reso- 
nances and depends on particle interactions - through the 
Luttinger parameter K - and on the impurity strength 
through a phase shift. At intermediate distances and 
for Luttinger parameter K > 1/2, the interaction is of 
the RKKY type, featuring Friedel oscillations with an 
envelope which is decaying as a power law with a renor- 
malized exponent. When 1/4 < K < 1/2, Friedel os- 
cillations are still present at finite T but are not decay- 
ing anymore. With ultracold fermionic atoms confined 
in a cigar-shaped trap, as described in detail in part I, it 
should be possible to detect this oscillating Casimir force. 



Ref. |l-y . 1 1 51 ] and are reviewed by Giamarchi [18j. In 
these works, the forward scattering term on the impu- 
rities was not taken into account. When this term is 
neglected, and in the absence of a gate voltage, the cri- 
terion for tunneling resonances in a Luttinger liquid is 
cos(k F r) =0 [3, \m, [H| . It is trivial to redo these calcu- 
lations starting from the action ^ including the phase 
shift 5 = —2Kg/u resulting from forward scattering on 
the impurities. The only modification amounts to replace 
2kp r by 2kpr + 8 and therefore the criterion for tunnel- 
ing resonances becomes coa(kpr + 8/2) — 0. As in the 
case of non- interacting fermions, this coincides with the 
position of the maxima in V\2 . 



Appendix B: DISAPPEARANCE OF THE 
OSCILLATIONS IN Vw 



Acknowledgments 

We acknowledge useful discussions with I. Safi, 
H. Saleur, N. Prokof'ev and P. Faccioli. We also thank 
V. Meden for very instructive discussions and comments 
on the present manuscript. 



Appendix A: TUNNELING RESONANCES 

We first consider a single fermion at the Fermi surface - 
with incoming wave vector kp - scattering on two delta 
impurities giS(x) and g25(x — r). Using the transport 
matrix formalism, it is easy to obtain its transmission 
probability across the double barrier structure 



P(k F ) = [1 
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+ 27172-y/ (1 + 7i 2 )(l + ll) cos(2k F r + 8) 



where jj = mgj/kp and 5 = — arctan7i— arctan72. This 
probability has maxima - called tunneling resonances - 
when cos(k F r + 5/2) = 0. On resonance, the energy of a 
quasi-stationary state in between the impurities matches 
the Fermi energy. In the case of symmetric impurities 
(71 = 72), there is perfect transmission on resonance. 
For asymmetric impurities, the transmission probability 
is always smaller than one. The condition for tunneling 
resonances coincides with the position of maxima of Vi 2 . 

We now turn to interacting fermions. Tunneling 
resonances in a Luttinger liquid were first studied in 



Here we discuss an interesting phenomenon in the 
interaction energy V12, which was pointed out to us 
by V. Meden [13]. Consider non-interacting spinless 
fermions on a lattice in a tight-binding model with hop- 
ping amplitude t, lattice spacing b, and filling factor / = 
pob = 1/2. The Fermi momentum is kp = nf/b = n/(2b) 
and the Fermi velocity is vp = 2tbsm(fir) — 2tb. This 
model is expected to have the same low energy behav- 
ior as the continuum model considered in the present 
paper. There are also two site impurities with dimen- 
sionless strength 71 and 72 separated by a distance r, 
which can only take discrete values r = integer x b. 
Therefore 2kpr = n /(2b) x integer x b = integer x n, 
which means that expi% = (— l) r / b exp id, where 5 — 
— arctan[(7i + 72)/(l — 7172)]- If in addition the impu- 
rities are fined tuned such that 7172 = 1, the phase shift 
5 = —ir/2 and Eq. ([TJ becomes: 
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This shows that the interaction appears to be purely at- 
tractive - as Li 2 (-l/4) < Li 2 (-(7i + 72)~ 2 ) < - and 
not oscillating, when probing inter-impurity distances 
r on the lattice only. This phenomenon occurs when 
three conditions are fulfilled: the system is half-filled 
(kp = 7r/(26)), the impurity separation only takes dis- 
crete values on a lattice (r = integer x b) and the impurity 
strengths satisfy 7172 = 1. For a numerical study of the 
half-filled lattice model of spinless fermions with nearest- 
neighbor interaction, and with two impurities, showing 
that this phenomenon survives in the case of interacting 
fermions, see Ref. [131 . 
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